Dynamical Criterion for Quantum Chaos : Entropy Production in Subsystems 
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We examine the conjecture that entropy production in subsystems of a given system can be used 
as a dynamical criterion for quantum chaos in the latter. Numerical results are presented for finite 
dimensional spin systems as also for the quantum baker's map. Of especial importance is the power 
spectrum of the entropy production which gets progressively more and more broad-banded as the 
degree of correlation in the Hamiltonian matrix is made to decrease. 
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More than twenty years after the discovery of dynami- 
cal localisation [1], possible signatures of quantum chaos 
are still under active exploration. While investigations 
have lent overwhelming support to a set of sharp con- 
jectures and predictions relating to stationary aspects of 
the problem like level statistics and eigenvector distribu- 
tions (see, e.g., [2-4]), corresponding dj/namica/ criteria of 
quantum chaos are only recently beginning to take shape. 
One recalls that a simple-minded search for chaotic fea- 
tures in quantum time evolution of a system along lines 
similar those in classical evolution should yield negative 
results in consequence of the linearity of the Schrodinger 
equation. A more fruitful approach, on the other hand, 
would be to focus on reduced states of subsystems. Re- 
cent investigations have focussed on the reduced entropy 
of subsystems [5] as an important indicator in the con- 
text of quantum chaos. For a set of relevant results on 
reduced entropy as a measure of correlation in subsys- 
tems, see [6,7] 

The stationary criteria of level statistics and eigenvec- 
tor distribution share the satisfactory feature that they 
are specifically quantal in nature, having no direct ref- 
erence to the semiclassical limit. Dynamical criteria like 
the spreading of wavepackets examined in earlier investi- 
gations (see, e.g., [8,9]) on quantum chaos are not entirely 
satisfactory on this count since they refer to time scales 
pertaining to the semiclassical limit. By contrast, the 
dynamic criterion of entropy production in subsystems 
retains its relevance in the deep quantum regime as well 
as in the semiclassical one. 

Quantum chaos has acquired a special relevance in the 
context of quantum information processing and quantum 
computation [10] in that it reflects in the manner in which 
decoherence and entanglement [11-13] arise in a system 
of qbits. Entropy production can be used in this con- 
text as a natural measure of decoherence, and it seems 
appropriate to look at various features of entropy produc- 
tion in order to formulate a sufficiently satisfactory and 
universal dynamical criterion for quantum chaos. The 



degree of entanglement, which is represented by entropy 
production under certain circumstances, has also been 
shown to be sensitive to the presence of classical chaos 
[6,5,14]. The present paper folows this line of approach 
and aims to point out that features of entropy produc- 
tion are relevant in characterising quantum chaos, one 
such useful feature being the time dependence of entropy 
production as revealed by its power spectrum. 

More precisely, starting from an initial state p(0) of a 
system S made up of subsystems 5*1 and S2, we let the 
system evolve for a time t and obtain the reduced state 
Pfl(i) of Si on taking partial trace over 82- We then look 
at the time dependence of pu. Even when one starts 
from a pure state p(0), pa turns out more often than not 
to correspond to a mixed state and there takes place, in 
general, a decoherence with the passage of time. As men- 
tioned above, one can measure the extent of decoherence 
by the reduced von Neumann entropy 



Sfl(i) = -trace{pR{t)lnpR{t)). 



(1) 



We compare features of time dependence of s_r(^) for sys- 
tems conforming to the predictions of random matrix the- 
ory (RMT) with the corresponding features for a 'regular' 
(see below) syestm. We present evidence that while the 
former is characterised by a broad-band power spectrum, 
the latter involves sparsely distributed peaks, indicative 
of a quasiperiodic time dependence. The system we con- 
sider is the simple one of an interacting assembly {S) of 
n spins or, in other words, an n-qbit system. The state 
space has dimension TV = 2", and the Hamiltonian de- 
termining the time evolution of S" is an iV x matrix. 
While quantum chaos in the true sense is expected to 
emerge only in the limit A^ — > 00, we consider finite (and 
rather low) values of A^ here so as to present in as simple 
a context as possible the features we aim to underline in 
this note. Finite dimensional considerations are, more- 
over relevant in numerous contexts relating to quantum 
computation. 
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Thus, wc consider, first, an N x N real symmetric 
Hamiltonian matrix with randomly distributed ma- 
trix elements for which the stationary features predicted 
by RMT arc conformed to and, next, a second Hamilto- 
nian matrix that is equivalent to the Harper Hamil- 
tonian on a torus. In the qbit representation the Hamil- 
tonians correspond to strong and long-range interactions 
among the qbits. Alternatively, either of these can also 
be interpreted to represent the Hamiltonian of a chain of 
N fermions of the form 

^aij{ci+Cj+CiCj+), (2) 

where = aji, and Ci, Ci^ denote fermion annihilation 
and creation operators. For He , the coefficients 
chosen as randomly distributed matrix elements (subject 
to reality ad hermiticity) while, for H^, we take 

Hr=llT + ^2V. (3) 

Here 71 , 72 are real constants, 
Ti,i+i = Ti+i,i = Tijv = Tni = 1/2 {i = 1, ...N - 1), 

(4) 

and 

Vij=cos{2TT3/N)6ij. (5) 

This represents a Harper's Hamiltonian 

Hr = 71 cos(27rp/P) + 72 cos{2wq/Q), (6) 

on a torus of area PQ with basic states {U = PQ /2'kN) 
and with torus boundary conditions Wi = 0^2 = 0. The 
parameters 71,72 have been set at (resp.) 0.5, 2.5 in our 
computations. In the fermion-chain interpretation, 
represents a Hamiltonian in the 1-fermion subspace for 
the chain [14]. However, in the following, we refer to 
He and H^ simply as Hamiltonians of an interacting spin 
system. 

With the Hamiltonian of the system S defined as 

above, we focus on a subsystem Si made up of, say, 
p(< n) specified spins and look at the reduced density 
matrx (RDM) given by 

PR{t)=Tr^^-^\e^p{0)e'^], (7) 

where H = He or Hr as the case may be, and where 
j'.j.l'i-p) iii([i(;.;ites partial trace with respect to states of 
the remaining {n—p) number of spins. We then calculate 
sii{t) of Eq.(l). We present below numerical results for 
n = S {N = 256), and p = 5, i.e., a subsystem made up of 
five spins. For convenience of presentation, we represent 
all the quantities as dimensionless; in other words, we use 
arbitrary but consistent units. Fig. 1 gives a comparison 
of the typical time-dependence of Sij(t) for (a) He and 
for (b) Hr with the same value of %. For the sake of com- 
parison we subtract a constant term from each matrix 
element of He, Hr so as to reduce to zero the mean of all 
the matrix elements. 
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FIG. 1. Entropy production as a function of time for (a) 
He, and (b) Hr, {h = .1 in both cases); note the erratic char- 
acter of the variation in (a) as opposed to the more regular 
variation in (b). 

Referring to Fig. 1 (see also Fig. 6) we note the fol- 
lowing distinctive features of Sii{t) : (i) there is an initial 
'transient' regime in which Sii{t) increases more rapidly 
for He (almost instantaneous in the scale of Fig. 1, but 
discernible in Fig. 6(a) for the quantum baker's map) 
than for Hr; (ii) after the initial transient regime there 
arises a 'steady' state in which Sij(t) fluctuates accord- 
ing to some invariant pattern; however, the mean value 
around which the fluctuation takes place is, in general, 
larger for He than for Hr', (iii) the fluctuation is more 
erratic for He than for Hr - while the latter appears to 
be quasiperiodic in nature, the former is random by com- 
parison. 




frequency frequency 

FIG. 2. Power spectrum of SR{t) for (a) He, and (b) Hr; 
the initial transient regime has been ignored while olrtain- 
ing the power spectrum; note the broad-band nature of the 
spectrum in (a) in contrast to (b). 

The contrast between the modes of fluctuation in the 
two situations under study is seen more clearly in Fig. 
2 where we present the power spectrum of the fluctua- 
tions in Sr. One finds that the power spectrum for He 
is a broad-band one while that for Hr has a relatively 
smaller number of sharp peaks (the power spectra have 
been obtained with the help of a standard FFT pack- 
age, and the symmetry about the mid-point of the fre- 
quency range has no physical relevance, being an artifact 



2 



resulting from the program). It may be mentioned that 
the quahtative features indicated above persist against 
changes in the value of U and the dimension of the sub- 
system considered. Details will be presented clsewhcr. 

In order to set in perspective the distinction between 
Hr and He one can compute the nearest neighbour level 
spacing distribution for the two Hamiltonians. One ob- 
serves in Fig. 3 that the level spacing distribution for He 
is close to the Wigner distribution implied by RMT, while 
that for Hr deviates to a large extent from the Wigner 
distribution and is, in fact, closer to the Poisson distri- 
bution followed by systems close to integrable classical 
systems. In other words, there exists a discernible corre- 
lation between stationary features implied by the RMT 
and dynamic features of entropy production in subsys- 
tems. Put differently, entropy production can indeed be 
taken as a valid dynamical criterion for quantum chaos - 
a criterion, moreover, that makes no direct reference to 
classical time evolution. While the latter feature is ap- 
parent from the manner in which the criterion has been 
formulated, one also finds support for this assertion from 
the fact that features of entropy production mentioned 
above appear even in situations with higher values of h 
(not shown in our figures) compared to the value used 
here. 




FIG. 3. Nearest neighbour level spacing distribution for (a) 
He, and (b) Hr; the dotted line compares the distribution in 
each case with a Wigner distribution with appropriate mean 
and variance . 

A quantum system represented by a Hamiltonian with 
randomly distributed matrix elements is, in a sense, anal- 
ogous to classical systems exhibiting 'hard' chaos. By 
contrast, analogs to 'regular' classical systems would be 
ones with Hamiltonian matrices having correlated ma- 
trix elements. In between, there exist Hamiltonian ma- 
trices with matrix elements having various intermediate 
degrees of correlation. It would be interesting to see if 
features of entropy production in subsystems of such sys- 
tems is, in some sense, intermediate between those of 
the systems considered above, analogous to the fact that 
systems characterised by 'soft' chaos in the semiclassi- 
cal limit exhibit level spacing distributions intermediate 
between the Wigner and Poisson distributions (see e.g., 
[15]). 



With this in view we modify our Hamiltonian Hr by 
replacing a given fraction of matrix elements with ran- 
domly distributed ones. To be more precise, we construct 
a family of Hamiltonian matrices H{f) depending on the 
parameter /(O < / < 1) where increasing values of / 
correspond to greater degrees of correlation between the 
matrix elenments. H{f) is obtained from Hr by replacing 
i-^r)mn with randomly distributed numbers for m, n sat- 
isfying |m — n| > fN, where / is some chosen fraction. 
Plots of sii{t) for decreasing values of /(not shown in our 
our figures) indicate that decreasing values of the degree 
of correlation indeed correspond to a gradual transition 
in features of entropy production from those of regular 
systems (Fig. 1(b)) to chaotic ones (Fig. 1(a)), being 
intermediate between the two. This is further corrobo- 
rated from Figure 4 where we display the corresponding 
power spectra for two chosen values of /. One observes 
that with decreasing / the power spectrum of entropy 
production becomes more and more broad-banded. 




frequency frequency 



FIG. 4. Power spectrum of SR{t) with (a) / = 0.9, and (b) 
/ = 0.8; in (a), the spectrum becomes more broad-banded 

with a decrease in correlation among the elements of the 
Hamiltonian matrix, corresponding to a decrease in /. 

The fact that the dynamic criterion of chaos indicated 
in this paper including, in particular, the power spec- 
trum of entropy production, is well correlated with the 
static criteria implied by RMT is borne out by looking 
at eigenvector distributions of H(f) for various values of 

The results are presented in Fig. 5 in which the resid- 
ual parameters [16] for the eigenvectors of H{f) are dis- 
played, once again for a chosen pair of values of /. The 
residual parameter for any given eigenvector is defined as 

r = ^ij:(P^-P'^^'W, (8) 

where the set of numbers pk represent the probability 
distribution for the elements of the eigenvector and pfe^"^ 
stands for a Gaussian distribution with the same mean 
and variance. Low values of r for all the eigenvectors im- 
plies a correspondingly high degree of allegiance to the 
predictions of the RMT. One observes from the figures 
that decreasing values of / indeed corresponds to an in- 
creasing degree of conformity to the stationary features 
implied by RMT. 
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FIG. 5. (a) Residual parameters of eigenvectors for / = 0.9 

(upper trace) compared with those for an uncorrelated Hamil- 
tonian (/ = 0, lower trace); note a large deviation from pre- 
dictions of RMT; (b) corresponding comparison for / = 0.7 
(upper trace) with the uncorrelated case (/ = 0); the two sets 
of residual parameters are now closer to each other. 



While entropy production in subsystems can be 
loooked upon as an 'internal' dynamical criterion for 
chaos, dynamical features of an 'external' nature have 
also been discussed in the Itcrature in various related con- 
texts [19-22], namely the distinctive features of chaotic 
systems in interaction with or perturbed by external 
ones. 

Other features of entropy production in subsystems, 
including fractal characteristics, will be presented in a 
future publication. 



The above features of reduced entropy for He are also 
in evidence for the quantum baker's map [17]. Fig. 6 
shows (a) the time-dependence of s^, and (b) the power 
spectrum, which are of the general nature of Fig.s 1(a), 
2(a) (as opposed to 1(b), 2(b)) respectively. Interesting 
details of entropy production for the quantum baker's 
map, as also the quantum cat map [18], both known to 
be chaotic in the classical limit, will be presented else- 
where. 
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FIG. 6. Entropy production in the quantum baker's map 
with dimension N = 128, the subsystem chosen being of di- 
mension 8; (a) reduced entropy as a function of t; note the 
fast rise in the transient regime; (b) power spectrum, showing 
an emerging broad-band structure for even the comparatively 
low value of N chosen . 

In summary, we have presented evidence that features 
of entropy production, including the power sprctrum in 
the steady state, in subsystems of a given system can be 
used as dynamical criteria of quantum chaos in the latter. 
Such dynamical criteria complement the stationary cri- 
teria discussed in the literature. The features presented 
show a progressive transition from regular to chaotic sys- 
tems when the degree of correlation in the underlying 
Hamiltonian matrix is made to decrease. 
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